Failure process of a bundle of plastic fibers 
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We present an extension of fiber bundle models considering that failed fibers still carry a fraction 
< a < 1 of their failure load. The value of a interpolates between the perfectly brittle failure 
(q = 0) and perfectly plastic behavior (a = 1) of fibers. We show that the finite load bearing capacity 
of broken fibers has a substantial effect on the failure process of the bundle. In the case of global 
load sharing it is found that for a -t 1 the macroscopic response of the bundle becomes perfectly 
plastic with a yield stress equal to the average fiber strength. On the microlevel, the size distribution 
of avalanches has a crossover from a power law of exponent w 2.5 to a faster exponential decay. For 
localized load sharing, computer simulations revealed a sharp transition at a well defined value a c 
from a phase where macroscopic failure occurs due to localization as a consequence of local stress 
enhancements, to another one where the disordered fiber strength dominates the damage process. 
Analysing the microstructure of damage, the transition proved to be analogous to percolation. At 
the critical point a c , the spanning cluster of damage is found to be compact with a fractal boundary. 
The distribution of bursts of fiber breakings shows a power law behaviour with a universal exponent 
w 1.5 equal to the mean field exponent of fiber bundles of critical strength distributions. The model 
can be relevant to understand the shear failure of glued interfaces where failed regions can still 
transmit load by remaining in contact. 

PACS numbers: 46.50.+a,62.20.Mk,64.60.-i 



I. INTRODUCTION 

The failure of heterogeneous materials under various 
types of external loading conditions has attracted contin- 
uous scientific and technological interest during the past 
decade [1, 2]. Both the macroscopic strength and the 
process of damaging of loaded specimens strongly depend 
on the disordered microscopic properties of the material. 
Hence, most of the theoretical studies are based on dis- 
crete models which can account for the disordered ma- 
terial properties and their interaction with the inhomo- 
geneous stress field naturally arising in a damaged speci- 
men. Fiber bundle models (FBM) are one of the most im- 
portant theoretical approaches in this field [3] , which also 
served as the basis for the development of more compli- 
cated micromechanical models of fracture [4-6] . In spite 
of their simplicity, FBMs capture the most important in- 
gredients of the failure process and make it also possible 
to obtain several characteristic quantities of high interest 
in closed analytic form. Based on FBMs important re- 
sults have been obtained for the macroscopic response of 
the loaded specimen [7], and for the temporal [8-11] and 
spatial structure of damage on the microlevel. In the 
framework of FBM the analogy of fracture and critical 
phenomena [12-15] has also been addressed, which is of 
high practical importance for the forecasting of imminent 
failure of loaded systems [9, 16]. 

Fiber bundle models have also been adopted to study 
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the failure of glued interfaces of solid blocks [17-20]. 
Such interfaces as a part of complex constructions, are 
assumed to sustain various types of external loads. In 
fiber reinforced composites, where fibers are embedded 
in a matrix material, the fabrication of the fiber-matrix 
interface strongly affects the mechanical performance of 
the composite. Very recently we have shown that under 
shear loading of glued solid blocks, the interface elements 
may suffer not only stretching but also bending [21]. In 
order to capture this effect we proposed to discretize the 
interface in terms of beams which can be elongated and 
bent, and break due to both deformation modes in a 
complex way. During the gradual failure of interfaces of 
solid blocks under shear, damaged regions of the inter- 
face can still transmit load contributing to the overall 
load bearing capacity of the interface. This can occur, 
for instance, when the two solids remain in contact at 
the failed regions and exert friction force on each other. 
In many applications the glue between the two interfaces 
has disordered properties but its failure characteristics is 
not perfectly brittle, the glue under shear may also yield 
carrying a constant load above the yield point. 

We present an extension of models of the shear fail- 
ure of glued interfaces considering that surface elements 
after failure still can have a certain load bearing capac- 
ity. The disordered interface is represented by a parallel 
set of fibers with random breaking thresholds and lin- 
early elastic behavior until failure. The broken fibers 
are assumed to carry a constant load which is a frac- 
tion < a < 1 of their failure load. Varying the value 
of a the model interpolates between the perfectly brit- 
tle (a = 0) and perfectly plastic (a — 1) constitutive 
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behavior of fibers. Based on analytic calculations and 
computer simulations, we show that the finite load bear- 
ing capacity of failed fibers has a substantial effect on 
both the macroscopic response and microscopic damage 
process of the fiber bundle. When the load redistribu- 
tion following fiber failure is short ranged, an interesting 
phase transition is revealed at a specific value of a. 



II. MODEL 

In order to model the shear failure of glued interfaces, 
we recently introduced a model [21] which represents the 
interface as an ensemble of parallel beams connecting the 
surfaces of two rigid blocks. The beams are assumed 
to have identical geometrical extensions (length I and 
width d) and linearly elastic behavior characterized by 
the Young modulus E. In order to capture the failure 
of the interface, the beams are assumed to break when 
their deformation exceeds a certain threshold value. Un- 
der shear loading of the interface, beams suffer stretch- 
ing and bending deformation resulting in two modes of 
breaking. The stretching and bending deformation of 
beams can be expressed in terms of a single variable, i.e. 
longitudinal strain e = Al/l , which enables us to map 
the interface model to the simpler fiber bundle models. 
The two breaking modes can be considered to be indepen- 
dent or combined in the form of a von Mises type break- 
ing criterion. The strength of beams is characterized by 
the two threshold values of stretching E\ and bending £r 2 
a beam can withstand. The breaking thresholds are as- 
sumed to be randomly distributed variables of the joint 
probability distribution p{e\, £2). The randomness of the 
breaking thresholds is supposed to represent the disorder 
of the interface material. After breaking of a beam the 
excess load has to be redistributed over the remaining in- 
tact elements. In Ref. [21] we presented a detailed study 
which demonstrated that the beam model of sheared in- 
terfaces with two breaking modes can be mapped into 
a simple fiber bundle model of a single breaking mode 
by an appropriate transformation of the fibers' strength 
disorder. 

In the present paper, we extend our model by taking 
into account that failed surface elements can still carry 
some external load increasing the load bearing capacity 
of the damaged interface. For simplicity, our study is 
restricted to discretize the interface in terms of fibers 
which could then be further generalized to beams [21]. 
A bundle of parallel fibers is considered with breaking 
thresholds <j th in the interval < a t h < a th ax wrtn 
a probability density p(a t h) and distribution function 
P{&th) — j^ th p((r' th )dcr' th . We assume that after the 
breaking of a fiber at the failure threshold a\ hl it may 
retain a fraction < a < 1 of its ultimate load <r\ h , i.e. 
it will continue to transfer a constant load cxa\ h between 
the surfaces. This assumption can be interpreted so that 
at the damaged areas of the interfaces the two solids still 
remain in contact exerting for instance a friction force 
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FIG. 1: Constitutive behavior of a single fiber: the fiber shows 
linearly elastic behavior up to the breaking threshold a\ h , then 
it keeps a fraction < a < 1 of the ultimate load aa\ h . 



which may contribute to the overall load bearing capac- 
ity. In many applications the glue between the two inter- 
faces has disordered properties but its failure characteris- 
tics is not perfectly brittle, the glue under shear may also 
yield. The constitutive behavior of single fibers is illus- 
trated in Fig. 1. Note that the load carried by the broken 
fibers is independent of the external load, furthermore, 
it is a random variable due to the randomness of the 
breaking thresholds. Varying the value of a, the model 
interpolates between the perfectly brittle failure (a = 0) 
and perfectly plastic (a = 1) behavior of fibers. The load 
stored by the failed fibers reduces the load increment re- 
distributed over the intact fibers, which strongly affects 
the process of gradual failure occurring under quasi-static 
loading of the interface. In the following we present a de- 
tailed study of the model system varying the strength of 
plasticity a. For the range of load sharing the two limit- 
ing cases of global and local load redistributions will be 
considered after failure events. 



III. TRANSITION TO PERFECT PLASTICITY 

Assuming global load sharing (GLS) after fiber break- 
ing, the constitutive equation of the interface can be cast 
into a closed form. At an externally imposed deformation 
e the interface is a mixture of intact and broken fibers, 
which both contribute to the load bearing capacity of the 
interface. Since the broken fibers retain a fraction a of 
their failure load, at the instant of fiber breaking only 
the reduced load (1 — a)a\ h is redistributed over the in- 
tact fibers. Since the fraction of fibers having breaking 
threshold in the interval [e, £ + cfe] can be obtained as 
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FIG. 2: Constitutive behavior 17(e) of the plastic fiber bundle 
for uniform (a, b, c) and Weibull distribution with m = 2 (d, 
e, /) at a = 0.2 (a, d), a = 0.5 (b, e) and a = 1.0 (c, /). The 
contribution of intact odfbm and failed fibers opi are also 
shown. Note that ctdfbm is identical with the constitutive 
curve of simple dry fiber bundle models. 



p(e)ds, the constitutive equation u(s) reads as 

a(e) = Ee(l - P(e)) + a [ Ee'p{e')de', (1) 

Jo 



CDFBM 



<7pi 



where the integration is performed over the entire load 
history. The first term labeled ctdfbm provides the load 
carried by the intact fibers, which corresponds to the 
classical dry fiber bundle (DFBM) behavior [3, 8, 22]. 
The constitutive law of DFBMs is recovered in the lim- 
iting case a — 0, when the complete load of the failed 
fiber is transferred to the remaining intact fibers of the 
bundle. In the second term api, which accounts for the 
load carried by the broken fibers, the integral is calcu- 
lated over the entire load history of the interface up to 
the macroscopic deformation s. It can be seen in Eq. (1) 
that the value of a controls the relative importance of 
the elastic and plastic terms influencing the macroscopic 
response <r(e) and also the microscopic damage process of 
the system. When a is increased, less load is transfered 
to the intact fiber and in the limiting case a — 1 failed 
fibers retain their entire load so no load transfer occurs. 
In this report, we explore the influence of the parameter 
a when it is tuned between these two extremal cases. In 
the following calculations the value of the fibers' Young 
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FIG. 3: (Color online) Simulations of stress controlled load- 
ing of a bundle of N — 1.6 x 10 6 fibers with Weibull dis- 
tributed breaking thresholds (A = l,m = 2). For clarity, the 
occurrence of macroscopic failure is indicated by vertical lines. 
Increasing a the constitutive behavior becomes perfectly plas- 
tic. 



modulus was set to unity E = 1. 

We note that the plastic fiber bundle model resembles 
up to some extent to the continuous damage fiber bun- 
dle model (CDFBM) worked out in Refs. [11, 23]. The 
main assumption of the CDFBM is that due to the ac- 
tivation of certain internal degrees of freedom, the fibers 
undergo a gradual softening process reducing their Young 
modulus in consecutive partial failure events. The fibers 
always remain linearly elastic but with a Young modulus 



E(k) 



'E, where the multiplication factor < a < 1 



describes the stiffness reduction in a single failure event 
and k denotes the number of failures occurred. If the 
fibers can fail only once (k = 1) and keep their stiffness 
value constant, the constitutive law of the system reads 
as 



a(e) = Ee{l - P(e)) + aEeP{e). 



(2) 



It was demonstrated in Refs. [11, 23] that increasing the 
number of times k the fibers can fail, the CDFBM devel- 
ops a plastic plateau, however, with a mechanism com- 
pletely different from the one considered here. 

It is instructive to consider two fundamentally different 
cases of disorder distributions -P(e), namely bounded and 
unbounded ones, where the largest breaking threshold 
crYh' x takes a finite value or goes to infinity, respectively. 
In this report, we focus on two specific realizations, i.e. 
a uniform distribution between and c"^ ax 
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0<a th < aT x , (3) 
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FIG. 4: Critical strain e c (a), and critical stress <r c (6) as a 
function of 1 — a for a Weibull distribution with A = 1, m — 2; 
symbols: simulation results, solid lines: analytical expressions 
Eqs. (8,9,10). 



and distributions of the Weibull type 

P(a th ) = l-e-^/» m , (4) 

are considered where A and to denote the characteris- 
tics strength and Weibull modulus of the distribution, 
respectively. For our study the Weibull distribution has 
the advantage that the amount of disorder in the failure 
thresholds can easily be controlled by the value of m. 

The functional form of the constitutive behavior cr(e) is 
shown in Fig. 2 for both disorder distributions Eqs. (3,4). 
It is interesting to note that for a < 1 there always exists 
a maximum of cr(e), just as in the case of DFBM. Under 
stress controlled loading conditions, macroscopic failure 
occurs at the maximum of a(e) so that the position and 
value of the maximum define the critical stress a c and 
strain e c of the bundle, respectively. It can be observed in 
Fig. 2 that the value of a c and e c are both higher than the 
corresponding values of DFBM indicating that the pres- 
ence of plastic fibers increases the macroscopic strength 
of the bundle. The decreasing part and the plateau of 
(7(e) can be realized under strain controlled loading con- 
ditions gradually increasing e. Under strain control the 
local load on the fibers is determined by the externally 
imposed deformation so that there is no load redistribu- 
tion after fiber failure. The fibers break one-by-one in the 
increasing order of their failure thresholds a\ h — Ee\ h . 
When the deformation e approaches the maximum value 
of the breaking thresholds e^ x = a 7 ^ h MX /E, all fibers 
must fail gradually so that the load of intact fibers ctdfbm 
tends to zero, while that of the broken fibers ap\ takes a 
finite asymptotic value 

<7P1 -> a = aE I e'p(e')de' = a (a th ) , (5) 
Jo 

where the integral is equal to the average fiber strength 
(ath)- When the strength of plasticity a is increased, the 
critical strain e c and stress a c , furthermore, the asymp- 
totic stress of the plateau a increase. The value of the 
critical deformation e c can be obtained by differentiating 
Eq. (1) with respect to e and calculating the root [22] 



from which the critical stress follows as a c — a(e c ). Eq. 
(6) implies that in the limiting case of a — ► 1 the critical 
strain e c tends to the maximum of the breaking thresh- 
olds £™£ x ', where P(e^ LX ) = 1. For the uniform distri- 
bution Eq. (3) we obtain 

e c = -^— hence, e c — ♦ 2e° c = eft™. (7) 
1 — a/Z a— >1 

Here e° denotes the critical strain of DFBM e° c = e™ h ax /2, 
which can be obtained by setting a = in Eq. (6). It fol- 
lows that for unbounded threshold distributions like the 
Weibull distribution, e c diverges so that perfect plastic- 
ity is only reached in the limit e c — > oo. The functional 
form of the divergence is not universal, due to the struc- 
ture of the third term on the left hand side of Eq. (6), 
e c depends on the specific form of p(e). For the Weibull 
distribution, e c as a function of a reads as 

/ /i \ 1 / m 

e c = e° c {l-ay 1/m , where e° c = A ( — J (8) 

for any Weibull exponent to. Parallel to this, the decreas- 
ing part and the plateau of the constitutive curve cr(e) 
disappear so that a c and a converge to the same finite 
value, which is the average fiber strength (a t h) 

a -> (a th ) and a c -> (<j th ) . (9) 

The average fiber strength (ath) can be determined as 

max i / 1 \ 

( atfc ) = £^_ and (a th ) = -T - (10) 
2 to y to / 

for the uniform and Weibull distributions, respectively. 
Here T denotes the Gamma function. 

In order to illustrate this behavior, Fig. 3 presents con- 
stitutive curves for Weibull distributed fiber strength ob- 
tained by computer simulations of stress controlled load- 
ing up to the critical point with A = 1 and to = 2. It is 
apparent that in the limiting case of a —> 1 the constitu- 
tive curve a(e) reaches a plateau, indicating a perfectly 
plastic macroscopic state of the system. The position of 
the maximum e c of the constitutive curves, i.e. the end- 
ing point of the curves, rapidly increases as a approaches 
1, while the value of the maximum a c tends to a finite 
value. In agreement with the analytic predictions Eq. 
(8), simulations confirmed that s c diverges as a power 
law whose exponent depends on the parameters of the 
strength distribution (see Fig. 4). 

Controlling the external stress, the constitutive curve 
of the system Fig. 3 can only be realized up to the max- 
imum, since at the critical load a c abrupt failure of the 
bundle occurs breaking all the surviving intact fibers in 
a large burst. The fraction <f> of fibers which break in 
the final burst causing global failure can be determined 
as (j> = 1 — P(e c (a)), which is illustrated in Fig. 5 as a 
function of 1 — a for the specific case of a Weibull distri- 
bution 



1 - P(e c ) - e cP (e c ) [1 - a] = 0, (6) 



0(a) = e -Vm(i-«). 



(11) 
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FIG. 5: Fraction of intact fibers 0=1 — P(s c (a)) at the point 
of macroscopic failure e c vs. 1 — a for a Weibull distribution 
of A = l,m = 2; circles: GLS simulation results, solid line: 
analytical solution Eq. (11). 
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FIG. 6: Analytic solution of the avalanche size distribution 
-D(A) at various different values of a. For a — > the usual 
power law distribution is recovered, whereas for a — > 1 an 
exponential decay of D(A) is obtained. For the specific cal- 
culations a Weibull distribution was used with m = 2. 



It can be observed that as the system approaches the 
state of perfect plasticity a — * 1, <f) tends to zero. This 
demonstrates that more and more fibers break before 
global failure occurs, and perfect plasticity is obtained 
when the strongest fiber fails at the maximum of cr(e) 
(compare also to Fig. 3). This argument also implies 
that for a — > 1, the difference of the microscopic damage 
process under stress and strain controlled loading dis- 
appears, the fibers break one-by-one without triggering 
avalanches of breakings. 



IV. AVALANCHES OF FIBER BREAKINGS 

Under stress controlled loading of the fiber bundle, the 
load dropped by a breaking fiber is redistributed over the 
intact ones. This load increment can give rise to further 
breakings which then may trigger an entire avalanche of 
failure events. The distribution D{A) of avalanche sizes 
A is an important quantity for the dynamical description 
of the loaded system. For the case of classical DFBMs 
(a = 0) under GLS conditions the avalanche size distri- 
bution D(A) can be obtained analytically [8, 24] as an 
integral, from which the asymptotic form of the distribu- 
tion for large avalanches proved to be a power law 



D(A) oc A 



-5/2 



OO. 



(12) 



The value of the exponent 5/2 is universal, it does not 
depend on the details of the disorder distribution of the 
failure thresholds [8, 24]. 

In order to obtain the analytical solution for the 
avalanche distribution in the presence of plastic fibers 



a / 0, we can follow the derivation of Rcfs. [8, 24], 
taking into account that the average number of fibers 
a{e 1 a)de which break as a consequence of the load incre- 
ment caused by a fiber breaking at the deformation e, is 
reduced by a factor of (1 — a) 



ep(e){l-a) 

a{s,a)de = — ^r-—de. 

1 - P(e) 



(13) 



Taking into account that the critical deformation e c 
where macroscopic failure occurs also depends on a, the 
avalanche size distribution D(A) can be cast in the form 



D(A) A 



A-l 



ec(a) 
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J a(e,a) A - 1 e- a(£ ' a)A x (14) 



[1 — a(e, a)]p(e) de. 

For the specific case of the Weibull distribution with an 
arbitrary modulus m the general equation Eq. (14) can 
be written in the form 

DIA - a) - (15) 

+ A? m(l - ale- 4 -] , 

where A c depends on the amount of disorder m and on 
the strength of plasticity a 



A c = A + 



1 



m(l — a) ' 



(16) 
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FIG. 7: Distribution D(A) of avalanches of size A for various 
values of a obtained by computer simulations for a system of 
N = 1.6 -10 7 fibers with Weibull distributed failure thresholds 
m — 2. Satisfactory agreement is obtained with the analytic 
results presented in Fig. 6. 



In Eq. (15) 7 denotes the incomplete Gamma function 
[41]. Two limiting cases can be distinguished in the solu- 
tion: first, for a — > the classical power law dependence 
Eq. (12) is recovered. This analytic solution is illustrated 
in Fig. 6 for a Weibull distribution with m = 2, where a 
power law of D(A) is apparent for a < 0.9. However, for 
the limiting case of a — > 1, we have to consider the behav- 
ior of the argument A c of the analytic solution Eq. (15). 
For a w 1 , there will be a regime of A values where the 
term l/(m(l — a)) dominates over A resulting in a faster 
decay of the distribution D(A) than any power. Still, for 
any values of a in the limiting case A >> A c (a), the 
usual mean field power law behavior Eq. (12) is asymp- 
totically recovered. Avalanche size distributions D(A) 
obtained from computer simulations at various different 
values of a are presented in Fig. 7. In a good quantitative 
agreement with the analytic predictions, the numerical 
results can be well fitted by a power law of exponent 5/2 
for moderate values of a. However, for a > 0.9 strong de- 
viations from the power law Eq. (12) can be observed for 
intermediate avalanche sizes 1 < A < 10 3 , which appears 
to be an exponential decay. Although in the analytical 
solution the asymptotic power law behavior is still visible 
for very large A, see Fig. 6, computer simulations in Fig. 
7 show solely a very steep decrease. It can be seen in 
the analytic solution in Fig. 6 that the relative frequency 
of avalanches of size A > O(10 3 ) is D = 0(1O~ 30 ) for 
a = 0.99, so it would require extremely large systems to 
count any such events. The size of the largest avalanche 
A max is plotted in Fig. 8 as a function of a. Obviously, 
A max is a monotonically decreasing function of a whose 



FIG. 8: Size of the largest avalanche in a GLS simulation with 
a Weibull distribution of m = 2. 



decrease gets faster in the regime where the distribution 
D(A) exhibits the crossover to the faster decaying form. 

An important consequence of the analytic solution Eqs. 
(15,16) is that the characteristic avalanche size where the 
crossover occurs from a power law to a faster decaying 
exponential form also depends on the amount of disorder, 
i.e. the stronger the disorder is, the larger the crossover 
size gets at a given a. 



V. LOCAL LOAD SHARING 

From experimental and theoretical point of view, it is 
very important to study the behavior of the plastic bun- 
dle when the interaction of fibers is localized. In the 
case of local load sharing (LLS) under stress controlled 
external loading conditions, the load dropped by the bro- 
ken fiber is redistributed in a local neighborhood of the 
fiber giving rise to high stress concentration in the vicin- 
ity of failed regions. Stress concentration leads to corre- 
lated growth of clusters of broken fibers (cracks), which 
plays a crucial role in the final breakdown of the sys- 
tem, i.e. macroscopic failure of the bundle occurs due 
to the instability of a broken cluster which then triggers 
an avalanche of failure events where all the remaining in- 
tact fibers break. This effect typically leads to a more 
brittle constitutive behavior of the system and the ap- 
pearance of non-trivial spatial and temporal correlations 
in the damage process [17, 23, 25, 26]. 

In the plastic bundle, after a fiber breaks it still re- 
tains a fraction a of its failure load <7 t h so that only 
the amount (1 — a)a t h is redistributed over the intact 
fibers in the neighborhood. It implies that the load bear- 
ing broken fibers reduce the stress concentration around 
failed regions giving rise to stabilization which also affects 
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FIG. 9: (Color online) Constitutive law a(e) of the LLS bun- 
dle obtained by computer simulations of a system of size 
L = 401 for several different values of a. The inset shows 
a magnified view of a(e) for the regime a < 0.4. For clarity, 
vertical lines indicate the location of macroscopic failure. For 
the breaking thresholds a Weibull distribution was used with 
m = 2. 




FIG. 10: The relative difference of the critical stresses 
and a^ LS of global and local load sharing systems as a func- 
tion of a for two different values of the Weibull modulus m. 
The vertical lines indicate the critical values of a, which were 
obtained in Sec. VC. 



the temporal and spatial evolution of damage during the 
loading process. 

In the following we consider a bundle of N fibers or- 
ganized on a square lattice of size L x L with periodic 
boundary conditions. The fibers are assumed to have 
Weibull distributed strength Eq. (4), where the value of 
A is always set to unity and for the Weibull modulus 
two different values are considered: m — 2 (large dis- 
order) and m — 4 (smaller disorder). After a failure 
event the load dropped by the broken fiber (1 — ot)a\ h is 
equally redistributed over the nearest and next-nearest 
intact neighbors in the square lattice, i.e. the local neigh- 
borhood of a broken fiber contains at most 8 intact sites. 
Stress controlled simulations have been carried out for 
system sizes ranging from L = 33 to L = 801 varying the 
strength of plasticity < a < 1. 



A. Macroscopic response 

It has been shown for DFBMs where broken fibers 
carry no load, that the macroscopic response of the bun- 
dle when the interaction of fibers is localized follows the 
constitutive law of the corresponding GLS system with 
a reduced critical strain and stress, i.e. the LLS bundle 
behaves macroscopically in a more brittle way than its 
GLS counterpart [23, 25, 26]. Figure 9 shows the consti- 
tutive curve of a plastic bundle of size L = 401 for several 
different values of a. It can be observed that for a « 



the constitutive curve exhibits the usual LLS behavior, 
i.e. the macroscopic failure is preceded by a relatively 
short non-linear regime and global failure occurs in an 
abrupt manner. The position of the macroscopic failure 
defines the value of the critical strain e^ LS and stress 
<J^ LS . It is very interesting to note that when a is in- 
creased, the LLS constitutive curves practically recover 
the behaviour of the corresponding GLS system, i.e. for 
a > 0.4 the macroscopic failure occurs when reaching the 
plateau of <r(e). 

The convergence of the LLS system to the GLS macro- 
scopic behavior is better seen in Fig. 10 where the relative 
difference of the critical stresses a^ LS (a) and a^ LS {a) of 
the global and local load sharing bundles is presented. It 
can be seen in the figure that there exists a threshold 
value a c of a above which the macroscopic response of 
the LLS bundle becomes very close to the corresponding 
GLS system, while below a c the constitutive behavior 
of the bundle changes continuously from the usual LLS 
response with a high degree of brittleness (a = 0) to 
the global load sharing behavior. It seems that at a c 
a continuous transition occurs between the two regimes. 
The transition indicates that as a consequence of the re- 
duction of stress concentration around failed fibers, the 
bundle can sustain higher external loads and is able to 
keep its integrity until the maximum of u{e) is reached. 
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FIG. 11: Avalanche size distributions D(A) obtained by com- 
puter simulations for the system size L = 401 with local load 
sharing, using Weibull distributed failure thresholds m = 2. 
The power law fit is demonstrated for a = 0.4. In the inset 
the largest avalanche A max is plotted versus a. 



B. Bursts of fiber breakings 




FIG. 12: Latest stable configuration in LLS simulations of a 
system of size L = 401, with a Weibull strength distribution 
m = 2 at different values of the control parameter a (a) 0.0, 
(b) 0.35, (c) 0.4, (d) 0.6. The strength of the largest cluster 
Poc in the lattices are (a) 0.003, (b) 0.097, (c) 0.517, (d) 
0.999. Broken and intact fibers are indicated by black and 
white, respectively. 



The evolution of the macroscopic response of the sys- 
tem with increasing a is accompanied by interesting 
changes of the damage process on the micro-level, char- 
acterized by the avalanches of fiber breakings and the 
cluster structure of failed regions. The avalanche statis- 
tics presented in Fig. 11 shows remarkable features. For 
a w 0, due to the high stress concentration around 
failed fibers, the LLS bundle can only tolerate small 
avalanches so that the avalanche size distribution D(A) 
decays rapidly. With increasing a the higher amount of 
load kept by broken fibers can stabilize the bundle even 
after larger bursts, hence, the cut-off of the distributions 
moves to higher values. It is interesting to note that also 
the functional form of the distribution D(A) changes, i.e. 
when a approaches a c the exponential cut-off disappears 
and the distribution becomes a power law 

£>(A) ~ (17) 

for large avalanches. The exponent [m of the power 
law was determined numerically as [i LLS = 1.5 ± 0.07, 
which is significantly lower than the mean field value 
^jGLS _ 2.5 [8]_ Increasing a above the critical point 
an exponential cut-off occurs and the power law regime 
of large avalanches gradually disappears. Comparing Fig. 
11 to the corresponding GLS results presented in Fig. 7, 
it is apparent that above a c the LLS distributions D(A) 
have the same functional form and follow the same ten- 
dency with increasing a as the mean field results. It can 



be concluded that the avalanche statistics presents the 
same transitional behavior between the local load shar- 
ing and mean field regimes as observed for the macro- 
scopic response. The same value of fi LLS was obtained 
numerically for m = 4, indicating the universality of the 
exponent with respect to the strength of disorder. The 
transition is more evident in the inset of Fig. 11, where 
the size of the largest avalanche A max is plotted as a 
function of a. The sharp peak indicates the transition 
point whose position defines a c , while in GLS the largest 
avalanche A max was a monotonically decreasing smooth 
function (compare to Fig. 8). 

C. Spatial structure of damage 

Gradually increasing the external load in the fiber bun- 
dle, the weakest fibers break first in an uncorrelated man- 
ner. Since the load is redistributed solely over the intact 
neighbors of the broken fiber, the chance of fiber break- 
ings increases in the vicinity of damage regions. This 
effect can result in correlated growth of clusters of bro- 
ken fibers with a high stress concentration around their 
boundaries. The larger the cluster is, the higher stress 
concentration arises. Global failure of the bundle occurs 
when, due to an external load increment, one of the clus- 
ters becomes unstable and grows until all fibers break. 
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FIG. 13: Distribution n of the size S of broken clusters in LLS 
simulations with a Weibull distribution m = 2, for different 
values of a. The spanning clusters were excluded from the 
distributions for a > 0.4. 



The spatial structure of the damage emerging when the 
interaction of fibers is localized can be characterized by 
studying the statistics and structure of clusters of broken 
fibers. Former studies of the limiting case of very local- 
ized interactions have revealed that the size of the largest 
cluster in the system is rather limited, furthermore, it 
is independent of the system size. Since the clusters are 
relatively small, merging of neighboring clusters does not 
occur frequently. The clusters themselves are found to be 
compact objects dispersed homogeneously over the cross 
section of the bundle [23, 27, 28]. 

In Fig. 12 the latest stable configuration of the bun- 
dle is presented just before catastrophic failure occurs 
at the critical load <r^ LS for several different values of 
a. For a « we note only small clusters of broken 
fibers as it is expected for LLS bundles (Fig. 12a). With 
increasing a, these clusters grow and adjacent clusters 
can even merge further increasing the typical cluster size 
(Fig. 12b). Around the critical value of a « 0.4, a span- 
ning cluster of broken fibers seems to appear (Fig. 12c), 
whereas for higher values of a > 0.4 almost all fibers have 
failed (Fig. 12d) already by the time the critical stress is 
reached. The existence of very large clusters is the di- 
rect consequence of the increased load bearing capacity 
of broken fibers. 

Clusters of broken fibers were identified in the square 
lattice using the Hoshcn-Kopclman algorithm. We eval- 
uated the distribution of cluster sizes n(S) in the last 
stable configuration just before macroscopic failure oc- 
curs. The behavior of n(S) shows again the transitional 
nature we have observed for other quantities. It can be 



seen again in Fig. 13 that a well defined a c exists which 
separates two regimes: for a < a c the clusters are small 
and n(S) has a steep decrease. Approaching a c , the clus- 
ter size distribution n(S) tends to a power law 



n(S) ~ S~ 



(18) 



where the value of the exponent was obtained as r = 
2.35 ± 0.08 which is higher than the corresponding expo- 
nent of 2d-percolation on a square lattice r = 187/91 « 
2.0549 [29]. Note that in the regime where spanning clus- 
ters exist (a > 0.4), the distribution n(S) contains only 
the finite clusters. 




a 



FIG. 14: Average cluster size S av = ni2/mi as a function of a 
for different system sizes L. Simulation results were obtained 
with a Weibull distribution m — 2. 



In order to characterize the evolution of the cluster 
structure when a is changed and to reveal the nature of 
the transition occurring at a c , we calculated the aver- 
age cluster size S av as the ratio of the second and first 
moments of the cluster size distribution 



m 2 
mi ' 



(19) 



The fc-th moment of the distribution n(S) is defined 
as 



m fe = ^S fe n(S)-Sl 



(20) 



where the largest cluster is excluded from the summa- 
tion. Figure 14 presents S function of a for dif- 
ferent system sizes ranging from L = 33 to L = 801. It 
can be seen in the figure that for each value of L the av- 
erage cluster size S av has a maximum at a well defined 
value of a, which becomes a sharp peak with increasing 
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FIG. 15: Finite size scaling of the average cluster size S av pre- 
sented in Fig. 14. The good quality collapse obtained enables 
us to determine the value of the critical exponents 7 and v 
with a relatively good precision. 



L, i.e. the peak becomes higher and narrower for larger 
systems. The observed behavior is typical for continuous 
phase transitions, where the position of the maximum 
defines the critical point of the finite size system. Based 
on the analogy to critical phenomena we tested the valid- 
ity of the scaling law S av ~ L 7 ^(j)((a — a c )L 1 / v ), where 
<j> denotes the scaling function of S av [29, 30]. The re- 
sults presented in Fig. 15 were obtained by varying the 
values of the critical point a c and of the critical expo- 
nent of the susceptibility 7, and correlation length v until 
the best data collapse was reached. It can be observed 
in Fig. 15 that in the vicinity of the critical point a c a 
good quality data collapse is obtained using the values 
a c = 0.385 ± 0.01, 7 = 2.0 ± 0.15, and v = 1.0 ± 0.1, 
where the critical exponents are only slightly different 
from the percolation exponents of 7 = 43/18 w 2.389 
and i/ = 4/3rj 1.33 in 2d [29]. 

At the critical point a spanning cluster of broken fibers 
occurs which is much larger than the other clusters. In 
order to characterize the strength of the spanning cluster 
we calculated the probability Px>(aO that a failed fiber 
belongs to the largest cluster. For percolation the quan- 
tity Pqo plays the role of the order parameter whose value 
distinguishes the phases of the system. Similarly to per- 
colation lattices, we find numerically a sharp rise from 
Poo = to Poo = 1 at a c « 0.4, see Fig. 16. When the 
system size L is increased P^ tends to a step function 
indicating that the transition becomes sharper. Assum- 
ing the scaling law Poo ~ L~@/ V ip((a — a c )L l l v ) of the 
order parameter for finite size systems, where ip denotes 
the scaling function and f3 is the order parameter expo- 




FIG. 16: Order parameter P x vs. a for several system sizes 
L with the Weibull index m = 2. 



nent [29, 30], we replotted the data in Fig. 17. The good 
quality of the data collapse was obtained with the pa- 
rameter values a c = 0.33 ± 0.01, (3 — 0.15 ± 0.06, and 
v = 0.95 ±0.1. Note that the value of v agrees well 
with the one determined by the finite size scaling of the 
average cluster size S av , larger deviations occur only for 
the critical point a c . The order parameter exponent [3 is 
compatible with the percolation value (3 = 5/36 « 0.13 



D. Random crack nucleation versus crack growth 

The failure mechanism of disordered materials and its 
relation to the amount of disorder has long been discussed 
in the literature [1, 2, 23, 28, 31-34]. When the mate- 
rial has a low degree of disorder only a small amount 
(if any) of damage occurs prior to macroscopic failure. 
In this case even the nucleation of a single microcrack 
can lead to localization and abrupt failure of the sys- 
tem. Increasing the amount of disorder, the macroscopic 
failure is preceded by a larger and larger precursory ac- 
tivity, i.e. a large amount of damage accumulates and 
local breakings can trigger bursts of breaking events [8] . 
Since cracks nucleate randomly, the process of damage 
before localization resembles percolation up to some ex- 
tent. Stress concentration builds up around failed regions 
which might lead to correlated growth of the nucleated 
cracks [12, 23, 31, 35]. Increasing the strength of dis- 
order, correlation effects become less dominating and in 
the limit of infinite disorder the damage accumulation 
process can be mapped to percolation [36]. 

We have shown above that in the plastic fiber bun- 
dle model (PFBM) , the load bearing capacity of broken 
fibers has a substantial effect on the process of failure 
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FIG. 17: Finite size scaling of the order parameter pre- 
sented in Fig. 16. The parameter values used to obtain the 
best quality data collapse agree well with the ones determined 
by the finite size scaling of S a v 



when the load redistribution is localized due to the reduc- 
tion of the stress concentration along cracks. In order to 
give a quantitative characterization of damage accumula- 
tion in our model, we determined the fraction of broken 
fibers pb at global failure a c as a function of the strength 
of plasticity a. The quantity pb can also be interpreted 
as the probability < pb < 1 that a randomly chosen 
fiber in the bundle is broken which makes it possible to 
compare the spatial structure of damage to percolation 
lattices [29] generated with the occupation probability 
p = Pb [1, 2, 26, 36]. The results are presented in Fig. 
18 for the system size L = 401 and Weibull parame- 
ters m = 2 and m = 4 plotting also the corresponding 
GLS results for comparison. In the case of local load 
sharing, when the failure load of fibers is almost entirely 
redistributed locally (a ~ 0) only a small damage can 
accumulate up to global failure Pb LS ~ 0.1 — 0.2 keep- 
ing the integrity of the system. Comparing the curves of 
different Weibull indices m it follows that the stronger 
the disorder is, the larger amount of damage the system 
can tolerate at the same value of a. In the vicinity of 
the respective a c , the breaking fraction p^ LS rapidly in- 
creases and converges to the maximum value p\ LS « 1, 
which implies that in the regime a > a c practically no 
localization occurs, the bundle can remain stable until 
almost all fibers break. 

It is instructive to compare this behavior to the case 
of GLS, where those fibers break up to the critical point 
whose breaking threshold falls below <r c , hence, pf LS (&) 
can simply be obtained as = P(a c (a)). It can be 

seen in Fig. 18 that under global load sharing for a ~ 
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FIG. 18: The fraction of broken fibers pt at a c as a function 
of a for fiber bundles of LLS and GLS with different strength 
of disorder m = 2 and m = 4. The vertical line indicates the 
critical point obtained as the position of the maximum of the 
average cluster size (see Fig. 14). The critical probability of 
percolation p c on the square lattice is indicated by the hori- 
zontal line. Note that for both disorder distributions in LLS, 
the location where p^ LS exceeds p c practically coincides with 
the corresponding critical point a c , indicating the percolation 
nature of the transition. 



a significantly larger fraction of fibers fails without de- 
stroying the system than in the LLS bundle. The break- 
ing fraction pb is a monotonically increasing function of 
a irrespective of the range of load sharing, however, in 
the vicinity of the critical point of LLS bundles ex- 
ceeds the smoothly rising GLS curves p^ LS ■ Note that 
depending on the threshold distribution P of fibers, even 
at a = the value of p^ LS can be smaller or larger than 
the critical percolation probability p c of the correspond- 
ing lattice type, since (contrary to fuse networks [31, 34] 
or discrete element models [37]) fracture in fiber bundles 
is not related to the appearance of a spanning cluster of 
failed elements. Varying a as a control parameter, for- 
mally the GLS results could be perfectly mapped onto a 
percolation problem: at the critical value of the control 
parameter af LS defined as P(a c (af LS )) = p c a spanning 
cluster occurs, which has a fractal structure, the average 
size of finite clusters has a maximum at the critical point 
and the cluster size distribution exhibits gap scaling [29] . 
However, this percolation is not related to the point of 
failure of the GLS bundle, the analogy to percolation is 
based purely on geometrical properties without any phys- 
ical relevance. 

Figure 18 shows that for localized load sharing the 
phase transition occurs when the damage fraction p^ LS 
reaches the critical percolation probability p c of the cor- 
responding lattice type. Due to the very localized load 
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sharing, only short range correlations arise in the system 
which arc further moderated by the finite load bearing 
capacity of broken fibers. Hence, in the vicinity of the 
transition point p^ LS (a c ) « p c holds and the evolution 
of the microstructure of damage shows strong analogy 
to percolation lattices. It can be seen in Table I that 
the critical exponents of the plastic fiber bundle model 
are slightly different from the corresponding exponents 
of percolation, furthermore, the usual scaling relations of 
percolation critical exponents [29] arc not fulfilled within 
the error bars. It has been shown for percolation that 
correlated occupation probabilities lead to the same criti- 
cal behavior as random percolation when the correlations 
are short ranged [38, 39], however, long range correlations 
result in changes of the critical exponents [38] . It is in- 
teresting to note that the value of the correlation length 
exponent v of PFBM is smaller than the value of ran- 
dom percolation which is consistent with the presence of 
relevant correlations [38]. We would like to emphasize 
that contrary to global load sharing, this percolation like 
transition has important physical consequences on the 
behavior of the fiber bundle. The failure process of the 
bundle is dominated by the competition of fiber break- 
ing by local stress enhancement due to load redistribu- 
tion and by local weakness due to disorder. Our detailed 
analysis revealed that the relative importance of the two 
effects is controlled by the parameter a. Below the crit- 
ical point a < a c high stress concentration can develop 
around cracks so that the failure of the bundle occurs 
due to localization. Above the critical point a > a c the 
macroscopic response of the LLS bundle becomes practi- 
cally identical with the GLS constitutive behavior show- 
ing the dominance of disorder. It is important to note 
that the size distribution of bursts of simultaneously fail- 
ing fibers becomes a power law at the critical point a c 
with an exponent [i equal to the value recently predicted 
for GLS bundles of so-called critical failure threshold dis- 
tributions [9, 10]. This can be explained such that the 
large avalanches of power law distribution occurring in 
the plastic fiber bundle model at a c (see Fig. 11) are dom- 
inated by the strong fibers of the bundle whose strength 
distribution is close to critical [9, 10]. 

The structure of the spanning cluster of the LLS bun- 
dle formed at the critical point a c has also remarkable 
features different from the spanning cluster of percolation 
[29]. The insets of Fig. 19 present representative exam- 
ples of the spanning cluster of a system of size L = 401 at 
two different disorder strengths. It can be observed that 
the clusters are compact, they practically do not have 
holes, there are no islands of unbroken fibers in the inte- 
rior of the cluster. This structure is a direct consequence 
of the merging of growing compact clusters where esp- 
cially large stress concentrations arise between the clus- 
ter surfaces breaking the fibers and filling the holes in 
the spanning cluster. We note that in the limiting case 
of very strong disorder a small amount of intact fibers 
may survive dispersed over the spanning cluster. The re- 
sult implies that the fractal dimension of the spanning 




l/lo 

FIG. 19: Structure of the spanning cluster at two different 
disorder strengths in a lattice of size L = 401. The perimeter 
length p(l) of the cluster is plotted as function of the length 
I of the yardstick normalized by the side length lo of the in- 
scribing square. The insets present the clusters analyzed. 



cluster of the LLS bundle is 2, which should be com- 
pared to the corresponding value of random percolation 
D = 91/48 w 1.896 where a finite amount of holes exists 
[29] even for short range correlated occupation probabili- 
ties [39]. The perimeter of the spanning cluster, however, 
has a fine structure, i.e. it has a large number of penin- 
sulas and valleys of all sizes. To reveal the structure of 
the perimeter, we measured its length p(l) as a function 
of the length of the yardstick I. It can be seen in Fig. 19 
that p(l) shows a power law dependence on I over almost 
two decades 

p(0~r 5 ', (21) 

where the value of the exponent proved to be S p — 
0.5 ± 0.03 for a Weibull distribution of fiber strength 
with m = 2. The power law Eq. (21) indicates that the 
perimeter line is a fractal with a dimension D p = l+5 p = 
1.5±0.03. The upper bound of the scaling range in Fig. 19 
can be attributed to the characteristic size of peninsulas 
of the spanning cluster, over which the rough structure 
of the perimeter disappears. Numerical calculations re- 
vealed that the fractal dimension of the cluster surface 
D p is not universal, i.e. it depends on the strength of 
disorder of the breaking thresholds. The insets of Fig. 19 
illustrate that a lower amount of disorder gives rise to a 
more regular, smoother cluster surface characterized by 
a lower value of D p . For the Weibull index m — 4 we 
obtained D p = 1.24 ± 0.05, which is significantly smaller 
than the corresponding value of m = 2. The surface of 
damage clusters should be compared to the hull of the 
spanning cluster of percolation with the fractal dimen- 
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TABLE I: Summary of the critical exponents of the plastic 
fiber bundle model with local load sharing. For comparison 
the value of the corresponding critical exponents of percola- 
tion are also shown. For the perimeter fractal dimension D p 
of PFBM a range is given. 



Critical exponents 


PFBM 


Percolation 


P 


0.15 ± 0.06 


5/36 w 0.13 


7 


2.0 ±0.15 


43/18 « 2.38 


r 


2.35 ± 0.08 


187/91 « 2.05 


V 


1.0 ±0.1 


4/3 « 1.33 


D 


2.0 


D = 91/48 « 1.896 


D p 


1.0 - 2.0 


7/4 = 1.75 


fi (Bursts) 


1.5 ±0.07 





sion D p = 7/4 = 1.75 [40] (see also Table I). 

VI. SUMMARY 

We introduced a fiber bundle model where failed fibers 
retain a fraction < a < 1 of their failure load. The 
value of the parameter a interpolates between the per- 
fectly rigid failure a — and perfect plasticity a = 1 of 
fibers. We carried out a detailed study of the effect of the 
finite load bearing capacity of fibers on the microscopic 
damage process and macroscopic response of fiber bun- 
dles considering both global and local load sharing for 
the load redistribution after fiber failure. Analytic cal- 
culations and computer simulations revealed that under 
global load sharing the macroscopic constitutive behav- 
ior of the interface shows a transition to perfect plasticity 
when a — ► 1, where the yield stress proved to be the av- 
erage fiber strength. Approaching the state of perfect 
plasticity the size distribution of bursts has a crossover 
from the mean field power law form of exponent 2.5 to a 
faster exponential decay. 

When the load sharing is localized it is found that 
the load carried by the broken fibers has a stabilizing 
effect on the bundle, i.e. it lowers the stress concentra- 
tion around clusters of failed fibers which has important 
consequences on the microscopic process of fracture and 



on the macroscopic response of the bundle. Extensive 
numerical calculations showed that at a specific value a c 
a very interesting transition occurs from a phase where 
macroscopic failure emerges due to stress enhancement 
around failed regions leading to localization, to another 
phase where the disordered fiber strength plays the dom- 
inating role in the damage process. 

On the macro-level, below the critical point a < a c 
the fiber bundle shows a brittle response, i. e. the macro- 
scopic failure is preceded by a weak non-linearity, while 
for a > a c the constitutive behavior of the LLS bundle 
becomes practically identical with the GLS counterpart. 
Analyzing the evolution of the micro-structure of damage 
with increasing a, the transition proved to be continuous 
analogous to percolation. Computer simulations revealed 
that the avalanche size distribution of fiber breakings be- 
comes a power law at the critical point with an universal 
exponent equal to the mean field exponent of bundles 
with critical strength distributions. The spanning clus- 
ter of failed fibers formed at the transition point proved 
to be compact with a fractal boundary whose dimension 
increases with the amount of disorder. The critical value 
a c is not universal, besides the lattice structure, it also 
depends on the strength of disorder. 

The plastic fiber bundle model can be relevant for the 
shear failure of interfaces where failed surface elements 
can remain in contact still transmitting load. Such glued 
interfaces of solids typically occur in fiber composites, 
where fibers are embedded in a matrix material. The fi- 
nite load bearing capacity of failed elements of the model 
can account for the frictional contact of debonded fiber- 
matrix interfaces and also for plastic behavior of the com- 
ponents. 
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